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CnI ■ Abstract. 

• I describe recent progress in terms of 3D hydrodynamical model at- 

^ [ mospheres and 3D line formation and their applications to stellar abun- 
dance analyses of late-type stars. Such 3D studies remove the free param- 

I eters inherent in classical ID investigations (mixing length parameters, 

■ macro- and microturbulence) yet are highly successful in reproducing a 
^ I large arsenal of observational constraints such as detailed line shapes and 

^ ■ asymmetries. Their potential for abundance analyses is illustrated by 

! discussing the derived oxygen abundances in the Sun and in metal-poor 

■ stars, where they seem to resolve long-standing problems as well as sig- 
I nificantly alter the inferred conclusions. 
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. 1. Introduction 

Q , Determining stellar element abundances play an integral role in most endeavours 

' to improve our understanding of stellar, galactic and cosmic evolution. The term 

c/2 , observed abundances is somewhat of a misnomer however, since the chemical 

' composition of course can not be inferred directly from an observed spectrum. 

J> . The obtained abundances are therefore never more trustworthy than the models 

employed to analyse the observations. The derivation of accurate abundances 
require realistic models of both the spectrum formation region and the spectrum 
I formation process, which is the topic of the present contribution. 

For a long time, understanding stellar convection has been one of the great- 
est challenges in stellar astrophysics. Convection can provide an efficient means 
of transporting energy and thereby greatly influence the stellar structure and 
evolution, while also mixing nuclear-processed material and preventing radiative 
diffusion and gravitational settling. For late- type stars, the surface convection 
zone reaches the stellar atmosphere, which directly affects the emergent spec- 
trum. The solar granulation is the observational manifestation of convection: 
concentrated, rapid, cold downdrafts in the midst of broad, slow, warm up- 
flows. Qualitatively similar granulation properties are expected in other solar- 
type stars, as indeed confirmed by 3D numerical simulations (e.g. Nordlund & 
Dravins 1990; Asplund et al. 1999; Asplund & Garcia Perez 2001; Allende Prieto 
et al. 2002) and indicated by observed line asymmetries. The up- and down- 
flows have radically different temperature structures (Stein & Nordlund 1998), 
which can not be approximated by normal theoretical ID hydrostatic model at- 
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mospheres with different effective temperatures T^^g (Fig. 1). In fact, neitfier 
tfie steep temperature structures of the ascending material nor the shallow gra- 
dients of the descending g ctS cirG well described by such ID models. Because of 
the photospheric inhomogeneities and the highly non-linear and non-local nature 
of spectrum formation, it is clear that no single ID model can be expected to 
properly describe all aspects of what is inherently a 3D phenomenon. 

Here I will describe recent progress in developing 3D hydrodynamical model 
atmospheres of late-type stars and some applications to stellar abundance analy- 
ses. In spite of their still inferior treatment of radiative transfer and line blanket- 
ing compared with existing ID model atmospheres based on the mixing length 
theory, detailed comparisons with a range of observational diagnostics suggest 
that these simulations are indeed highly realistic. 

2. 3D Hydrodynamical Model Atmospheres 

The 3D model atmospheres which form the basis of the abundance analyses 
presented here have been computed with a 3D, time-dependent, compressible, 
explicit, radiative-hydrodynamics code developed to study solar and stellar sur- 
face convection (Stein &; Nordlund 1998). The hydrodynamical equations for 
conservation of mass, momentum and energy are solved on a non-staggered Eu- 
lerian mesh with gridsizcs of ~ 100'^. The physical dimensions of the grids are 
sufficiently large to cover many (> 10) granules simultaneously in the horizon- 
tal direction and about 13 pressure scale-heights in the vertical. In terms of 
continuum optical depth the simulations extend at least up to log Tj^oss ^ — ^ 
which for most purposes are sufficient to avoid numerical artifacts of the open 
upper boundary on spectral line formation. The lower boundary is located at 
large depths to ensure that the inflowing gas is isentropic and featureless, while 
periodic horizontal boundary conditions are employed. The temporal evolution 
of the simulations cover several convective turn-over time-scales to allow ther- 
mal relaxation to be established and to obtain statistically significant average 
atmospheric structures and spectral line profiles. No magnetic fields have been 
included in the simulations described here although the code is capable of sim- 
ulating magneto-convection. 

In order to obtain a realistic atmospheric structure, it is crucial to have 
the best possible input physics, and properly account for the energy exchange 
between the radiation field and the gas. The adopted equation-of-state is that of 
Mihalas et al. (1988), which includes the effects of ionization, excitation and dis- 
sociation. The continuous opacities come from the Uppsala package (Gustafsson 
et al. 1975 and subsequent updates) while the line opacities are from Kurucz 
(1998, private communication). The 3D radiative transfer is solved at each time- 
step using 9 inclined rays (2 /x-angles and 4 (^-angles, plus the vertical) under the 
assumptions of local thermodynamic equilibrium (LTE, S\ = B\) and opacity 
binning (Nordlund 1982). The four opacity bins are designed to correspond to 
continuum and weak, intermediate and strong lines. The assignment of the orig- 
inal 2748 wavelength points into the different opacity bins follows from detailed 
monochromatic radiative transfer calculations of the ID averaged atmospheric 
structure. The opacity binning thus includes the effects of line-blanketing in a 
manner reminiscent of opacity distribution functions. In order to improve the 
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Figure 1. Temperature structure of a 3D solar model atmosphere 
illustrating the vastly different temperature gradients in the up- and 
downfiows. Vertical columns are connected with solid lines. 

numerical accuracy, the radiative transfer is solved on a finer grid and the re- 
sults subsequently interpolated back to the original grid for the hydrodynamical 
variables at each timestep. 

It is important to realise that the simulations contain no free parameters 
which are tuned to improve the agreement with observations. The adoption of the 
numerical and physical dimensions of the simulation box is determined by prac- 
tical computional time considerations, the need to resolve the most important 
spatial scales and the wish to place the artificial boundaries as far as possible 
from the region of interest. The code is stabilized using a hyper- viscosity dif- 
fusion algorithm with parameters specified from standard hydrodynamical test 
cases. These parameters are not changed in the actual stellar convection simu- 
lation and are therefore not freely adjustable parameters. It has been verified 
that the resulting atmospheric structures are insensitive to the adopted effective 
viscosity at the current highest affordable numerical resolution (Asplund et al. 
2000a). The input parameters discriminating different models are the surface 
gravity logg, metallicity [Fe/H] and the entropy of the inflowing material at the 
bottom boundary. The effective temperature of the simulation is therefore a 
property which depends on the entropy structure and evolves with time around 
its mean value following changes in the granulation pattern. 

Further details of the 3D hydrodynamical model atmospheres are available 
in Stein k Nordlund (1998), Asplund et al. (1999, 2000a,b) and Asplund & 
Garcia Perez (2001). 
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Figure 2. Predicted spatially resolved profiles of the Fel 608.2 nm 
line in the solar simulation. Note the bias of stronger lines in the 
upfiows due to their steeper temperature gradients in general. 



3. Spectral Line Formation 
3.1. General Description 

The 3D hydrodynamical model atmospheres described in the previous section 
form the basis for the 3D spectral line formation calculations. The Doppler shifts 
introduced by the convective and oscillatory motions cause line asymmetries and 
shifts as well as broaden the lines significantly. Spatially resolved profiles thus 
come in an amazing assortment as a result of the atmospheric inhomogeneities 
and velocities together with their correlations (Fig. 2). The concepts of micro- 
and macroturbulencc, which arc introduced in ID analyses in order to account 
for the missing line broadening, are not necessary in 3D calculations, as evident 
from the excellent agreement between predicted and observed line profiles (Fig. 
3, see also Asplund et al. 2000b) and the lack of significant trends in derived 
abundances with line strength (Asplund et al. 2000c; Asplund 2000). The main 
component of microturbulence does not arise from microscopic turbulent mo- 
tions but from gradients in convective motions that are resolved with the present 
simulations, which is also the main physical explanation for macroturbulence. 

The original 3D hydrodynamical simulations which includes deep, extremely 
optically thick layers are interpolated to a finer depth scale which extends only 
down to log tross ~ 2.5 prior to the spectral line calculations for improved numer- 
ical accuracy. At the same time, the simulations are interpolated to a coarser 
horizontal grid to ease the computational burden. Extensive tests have en- 
sured that this practice is acceptable even when considering fine details such 
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as line asymmetries; for abundance analysis purposes the procedure introduces 
non-noticable differences (< 0.01 dex for individual snapshots and much less for 
temporal averages). Depending on whether the 3D line formation is computed 
in LTE or in non-LTE, different assumptions and approximations are necessary, 
which are described next. 

3.2. 3D LTE Line Formation 

With the assumption of LTE {Sx = Bx), the level population is uniquely deter- 
mined by the local gas temperature from the Boltzmann and Saha distributions. 
With the source function known, it is then straightforward to solve the radia- 
tive transfer equation, which typically is performed for about 100 wavelength 
points across the line profile for 10-20 inclined outgoing rays for more than 50 
different snapshots. This allows asymmetries and shifts to be determined to 
an accuracy on the ms~^-level. For abundance determinations even individual 
snapshots give results accurate to within about 0.03 dex, since at any instant 
many granules are covered by the simulation box. All in all, a single tempo- 
rally and spatially averaged flux profile in 3D correspond most of the time to 
-^t * -^x * -^y * ^angles * -^A ^ 10^ ID radiative transfer calculations. In addition, 
each 3D profile is normally computed for at least three different abundances to 
enable interpolation to the requested line strength. Even then, such 3D LTE line 
calculations are achievable on current workstations thanks to efficient numerical 
algorithms and implementations. 

No micro- and macroturbulence enter the calculations but for flux profiles 
the projected stellar rotational velocity Wrotsini causes line broadening in ad- 
dition to the natural, pressure and convective broadening. Flux profiles are 
obtained after the actual radiative transfer calculations from a disk-integration 
of the stored profiles of the individual inclined rays for a given frotsini (Dravins 
& Nordlund 1990). The detailed comparison with observations can then pro- 
ceed through a x^-analysis by interpolating the computed profiles with varying 
abundances and adopting different Vrotsini. 

3.3. 3D Non-LTE Line Formation 

Recently, even detailed 3D non-LTE line formation calculations have become 

feasible (Kiselman 1997; Uitenbroek 1998; Botnen & Carlsson 1999; Asplund et 
al. 2002). Without the much simplifying assumptions of Saha and Boltzmann 
distribution for the level population, the rate equations assuming statistical equi- 
librium (dni/dt = 0) must be solved simultaneously with the radiative transfer 
equation. Due to the much increased computational effort in achieving this, 
sacrifies must be made in other areas, most notably in terms of number of snap- 
shots the calculations can be made for. With individual snapshots it is still 
feasible to carry out calculations for atoms consisting of hundreds of radiative 
transitions, each with 50-100 wavelength points. Fortunately, with 3D model 
atmospheres the non-LTE effects arc expected to be quite similar between dif- 
ferent snapshots due the relatively large area coverage always containing many 
up- and downflows, which is also confirmed by detailed calculations. The by far 
dominating part of 3D non-LTE computations of reasonably sized model atoms 
is the actual radiative transfer solution from the current estimate of the source 
function. This formal solution of the radiative transfer equation is performed 
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Figure 3. Upper panel: Comparing the predicted temporally and 
spatially averaged intensity profile of the Fe i 608.2 nm line (diamonds) 
with the observed solar profile (solid line). Lower panel: Corresponding 
ID profile using the Holweger-Miillermodel atmosphere with an opti- 
mized macroturbulence (Gaussian since radial-tangential macroturbu- 
lence is not applicable for intensity profiles) and a microturbulence of 
l.Okms"^. Note the inability of the ID case to predict the line shift and 
asymmetry, as well as the different required Fe abundance (+0.17 dex) 
for the same line strength. 
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for all wavelengths of all transitions for typically 24-48 inclined rays to compute 
mean intensities and flux profiles. The final converged solution is then obtained 
by a iterative procedure, normally employing an accelerated lambda iteration 
technique or similar, allowing convergence typically within 10 iterations. For 
small model atoms, the calculations are possible to perform on modern worksta- 
tions but for more extended atoms the use of supercomputers are better suited, 
not the least due to vast memory requirements with 3D non-LTE codes such as 
MULTI3D (Botnen & Carlsson 1999; Asplund et al. 2002). 

4. Stellar Abundance Analyses with 3D Model Atmospheres 

The realistic 3D hydro dynamical model atmospheres of late-type stars together 
with corresponding LTE and non-LTE line formation tools described above have 
a natural application to stellar abundance analyses. For abundance studies, the 
different mean photospheric structure between 3D and ID models, the existence 
of atmospheric inhomogenities and the broadening of spectral lines due to con- 
vective motions are all important ingredients for the line formation. To illustrate 
the possible impact of the new generation of 3D model atmospheres on abun- 
dance analyses, I will describe the situation for the most abundant cosmic metal, 
oxygen, in two important types of stars: the Sun and metal-poor halo stars. In 
both cases it seems like the adoption of 3D model atmospheres drastically affects 
the conclusions. 

4.1. Oxygen in the Sun 

The solar oxygen abundance ought to be well determined but astonishingly 
the commonly adopted value may be in error by almost a factor of two. For 
the Sun, five different types of diagnostics are potentially available: permitted 
high-excitation O I lines, the forbidden [O l] lines at 630 nm, the OH vibrational- 
rotational and pure rotational lines in the IR, and finally the OH electronic 
transitions in the UV (which have not been included in the present analysis). 
All have quite different sensitivities to the atmospheric structure and the line 
broadening. In addition, all have their pros and cons: the [Oi] lines are very 
weak and likely blended (Allende Prieto et al. 2001); the Ol lines are sensitive 
to departures from LTE (Kiselman 1993); the OH IR lines arc, like all molecular 
features, very temperature sensitive (Asplund & Garcia Perez 2001); the OH 
UV lines may be affected by missing opacity (Balachandran & Bell 1998) and 
possibly non-LTE effects (Asplund & Garcia Perez 2001). Unfortunately the 
different lines imply wildly discortant abundances, the difference amounting to 
about 0.3 dex both for the Holweger-Miiller (1974) and MARCS (Asplund et al. 
1997) ID model atmospheres. Normally the conclusions from the OH IR have 
been embraced (logeo ~ 8.90), largely due to the remarkably small line-to- line 
scatter (Sauval et al. 1984). An application of our new 3D model atmospheres 
to the problem, however, suggests that the solar oxygen abundance may be 
significantly lower (logeo ~ 8.65) than previously thought (Allende Prieto et al. 
2001; Asplund et al., in preparation). 

[O l] : Due to the exceptionally good agreement between predicted and ob- 
served line profiles in general with 3D model atmospheres (Fig. 3), disturbing 
blends are readily identified. This statement naturally presumes that detailed 
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comparisons have already been performed for a range of lines with similar prop- 
erties and typical line formation regions, which fortunately is the case for [Ol]. 
The very poor agreement between the theoretical and the observed profile, im- 
mediately signals that the [Oi] 630.0 nm line is indeed blended, which have 
occasionally been suspected but mostly thought to be relatively unimportant. 
The effect of the blending Ni I line reduces the derived abundance significantly 
(— 0.13dex), while the transition from ID to 3D model atmospheres and a re- 
vised (//-value arc of less importance (—0.08 and — 0.03dcx, respectively) but 
compounds the overall difference. The forbidden line implies a solar oxygen 
abundance of logeo = 8.69 ± 0.05 (Allcndc Pricto ct al. 2001). 

O I: The O I triplet around 777 nm is well-known to be susceptible to signif- 
icant departures from LTE. We have therefore investigated its line formation in 
detail with the 3D non-LTE code multiSd (Botnen k Carlsson 1998; Asplund 
et al. 2002), employing two independent snapshots from the 3D solar simulation 
and a model atom consisting of 23 levels and 65 radiative transitions (Kiselman 
1993). The adopted model atom is sufhciently comprehensive for this particular 
problem since even a two-level approach catches the essentials of the line forma- 
tion. Collisional excitation and ionization due to electrons are accounted for but 
not the corresponding case of hydrogen collisions, since the available classical 
recipes (Drawin 1968) likely over-estimates their importance by several orders of 
magnitude. Our 3D non-LTE calculations indeed reveal pronounced departures 
from LTE, in terms of abundance about 0.2 dex. The dominant non-LTE mech- 
anism is the same as with ID models, namely photon losses in the line itself. 
Preliminary the triplet implies a solar oxygen abundance of log eo = 8.67 it 0.03 
(Asplund et al., in preparation). 

OH vibrational and rotational lines: The slightly cooler temperature 
structure in our 3D solar model compared with the standard Holweger-Miiller 
(1974), which is normally the ID model of choice for solar abundance determi- 
nations, directly translates to a lower derived oxygen abTindancc from the OH 
IR lines. Equally important, however, is the presence of temperature inhomo- 
geneities due to the non-linear line formation. It therefore comes as no sur- 
prise that ID analyses typically will significantly over-estimate the abundances 
derived from molecular lines. We have reanalysed 70 and 127 OH vibration- 
rotation and pure rotation lines, respectively, which have previously been stud- 
ied by Grevesse et al. (1984) and Sauval et al. (1984) using up-to-date molecular 
line data. The 3D line formation calculations are based on the assumption of 
LTE, which is likely a reasonable approximation in view of the efficient colli- 
sional thermalization of the relevant levels. A potential source of error comes 
from the assumed chemical equilibrium in the OH molecule formation, which 
would tend to artificially decrease the oxygen abundance slightly. The magni- 
tude of this effect is not known at this stage but will be investigated further 
in the near future. Keeping this in mind, preliminary unweighted means yield 
logeo = 8.61 ± 0.03 for the vibrational lines and logeo = 8.67 ± 0.03 for the 
rotational lines. 

Summary: Although the exact value for the solar oxygen abundance de- 
pends on the relative weights the different diagnostics are given, it is clear that 
our detailed 3D investigation unanimously points to a significantly lower solar 
oxygen abundance than previously thought based on classical ID model atmo- 
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Figure 4. Upper panel: Derived O abundances from the OH 

vibrational-rotational lines in the IR using a 3D hydrodynamical sim- 
ulation of the solar atmosphere. The solid line denotes the least- 
square-fit to the data. Lower panel: Corresponding results for the ID 
Holwcgcr-Miillcr scmi-cmpirical model atmosphere. Note the distinctly 
higher derived abundances compared with the 3D case. 
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spheres. Preliminary we estimate logeo = 8.65 ± 0.05 which is lower by almost 
a factor of two. Such a radical re-evaluation for the most abundant metal in 
the star we know the best is certainly a sobering realisation when discussing 
derived stellar abundances based on inferior spectra and where the modelling 
uncertainties can be much larger. Our revision directly impacts also on stellar 
evolution calculations and helioseismology among other topics. 

4.2. Oxygen in Metal-Poor Stars 

The oxygen abundances in stars born at different epochs and locations carry 
vital information about the formation and evolution of galaxies. In particular 
the temporal (with metallicity as a practical although not perfect proxy) evolu- 
tion of the oxygen-to-iron ratio can yield estimates of the time-scales involved 
(O mainly from SNell while Fe largely from SNela) as well as give important 
clues to the nature of the elusive Pop III stars. It has long been known that halo 
stars show an over-abundance of oxygen ([0/Fe] > 0) but an exact estimate has 
been notoriously difficult to achieve due to the large discrepancies between the 
different abundance indicators. Typically the [O l] line implies a more or less 
plateau at [0/Fc] « +0.4 (Barbuy 1988; Nisscn ct al. 2002) but due to weak- 
ness of the line its use for [Fe/H] < —2.0 is restricted to analysis of giants and 
subgiants. Recent analyses of the OH UV lines in halo dwarfs on the other hand 
seem to suggest an almost linear trend with metallicity, reaching [0/Fe]~ +1.0 
at [Fe/H]= -3.0 (Israelian et al. 1998; Boesgaard et al. 1999) while the OH IR 
vibrational lines in giants again more support the conventional plateau (Melen- 
dez & Barbuy 2002). The confusing situation is completed by the fact that the 
oxygen abundances derived from the high-excitation O i triplet at 777 nm tend 
to lay in between the OH UV and the [O l] results. 

Wc have recently applied the new generation of 3D hydrodynamical model 
atmospheres of metal-poor stars (Asplund et al. 1999) in an attempt to investi- 
gate whether a possible resolution to the oxygen conundrum can be found. Due 
to the lack of radiative heating from spectral lines in halo stars, the adiabatic 
cooling from the rising gas in the upflows dominate the energy balance. Conse- 
quently, the temperature is much below the radiative equilibrium value in the 
optically thin layers, which is enforced in ID models. 

OH UV: As noted before, molecule formation is extremely temperature 
sensitive. The lower temperatures in 3D models therefore directly translates 
to much larger number densities of molecules and thus to stronger lines. The 
application of 3D model atmospheres appears to bring down the high advocated 
linear trend to something which better resembles the canonical plateau (Asplund 
& Garcia Perez 2001), revealing a potential and very significant systematic error 
in previous standard analyses using ID model atmospheres. This 3D LTE cor- 
rection can easily amount to 0.6 dex or more at [Fe/H] = —3.0. The final word 
on this topic is not written yet however, as the 3D OH analyses sofar make the 
crucial assumption of LTE both in the line formation and in the OH molecule 
formation. Existing non-LTE calculations for molecules arc very rare. Asplund 
& Garcia Perez (2001) found the possibility of significant departures from LTE 
in the OH line formation but with little metallicity dependence using a very 
simplified two- level approach in ID; the 3D case still remains to be investigated. 
In addition, due to the rapid cooling of the upflowing material it is not unlikely 
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that the time available is not sufficient for equilibrium chemistry to be estab- 
lished, in particular for metal-poor stars with their very low temperatures in the 
upper layers. 

[Oi]: One could imagine that the [Oi] line should be insensitive to the 
low surface temperatures in 3D models of halo stars insofar it originates from 
the ground level of the ionization stage where essentially all oxygen is residing. 
In reality though there are significant downward 3D corrections also for this 
line, since the line strength depends on the ratio of line and continuous opacity. 
The latter is dominated by H~ whose number density is proportional to the 
free electrons which are much reduced in the low temperature gas (Nissen et al. 
2002). In terms of abundance the 3D corrections amount to about — 0.2dex at 
[Fe/H] ~ —2.0 for dwarfs and subgiants. Whether a similar revision is needed 
for giants have to await the construction of 3D models for such stars but slightly 
smaller corrections would not be unexpected. 

O I: The most problematic of the available diagnostics in metal-poor stars 
is likely the O I lines. Due to their high excitation potential, the lines are formed 
in very deep layers where the differences between 3D and ID models are small. 
Already in ID analyses, the triplet yields higher O abundances than [O i] which 
will be aggrevated unless there are large metallicity-dependent 3D non-LTE 
effects on the triplet or the adopted Tcff-scale needs significant revision. Pre- 
liminary 3D non-LTE calculations for the halo subgiant HD 140283 (Asplund 
et al., in preparation) suggest quite similar abundance corrections to the solar 
case, ~ — 0.2dex, which is apparently not enough (Nissen et al. 2002). The 
best Teff-scale for halo stars likely comes from the IR flux method, which gives 
too low Teff to resolve the discrepancy even with 3D models (Asplund Sz Garcia 
Perez 2001). The use of 3D models may have resolved the problem between the 
[Ol] and the OH UV lines but clearly not all answers are in for the Ol triplet. 

Summary: It seems like the adoption of the new generation of 3D hy- 
drodynamical model atmospheres holds a crucial key to the resolution of the 
long-standing problem of oxygen abundances in halo stars. Although the final 
answer has not yet emerged, the existing results together with the recent down- 
ward revision of the solar O abundance described above point to a nearly flat 
[0/Fe] trend but with a slow rise for [Fe/H] ^ —2, i.e. somewhere in between 
the hotly contested plateau and linear trends (Nissen et al. 2002). 

5. Concluding Remarks 

The above examples illustrate that existing ID analyses may well be hampered 
by significant systematic errors which can severely impact the inferred conclu- 
sions. With the advent of a new generation of telescopes and instrumentation 
it is clear that the dominant source of error is no longer of observational nature 
but rather stems from uncertainties inherent to the actual analyses. There is 
obviously a great need for investing also significant resources in improving the 
modelling aspects in order to place the findings on firmer footing. It is notewor- 
thy that the main obstacle for achieving this is simply a lack of manpower, in 
spite of the existence of a very large community of astronomers relying on such 
model atmospheres for their applications. 
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